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The OA ansatz has attracted much attention recently, infinite-dimensional Kuramoto model could 
collapses to a two-dimensional system of order differential equations with it. In this paper, we 
propose the ensemble order parameter (EOP) equations to describe the dynamics for networks with 
a finite size. To verify the effectiveness of this method, we apply it into the star network and star- 
connected network. In the star network, numerous phase transitions among different synchronous 
states are observed, three processes of synchronization, one process of de-synchronization and a group 
of hybrid phase transitions, the processes of those transitions are revealed by the EOP dynamics 
and other nolinear tools such as time reversibility analysis and linear stability analysis. Also in the 
star-connected network, the two-step synchronization transition is observed. The process of it is 
still be revealed by the similar methods in the single star network. 


I. INTRODUCTION 

Understanding the intrinsic microscopic mechanism 
embedding in collective macroscopic behaviors of pop¬ 
ulations of coupled units in heterogenous network has 
become a focus in a variety of fields, such an biological 
neurons circadian rhythm, chemical reacting cells, and 
even society systems m- Numerous different emerging 
macroscopic states/phases have been revealed, and vari¬ 
ous non-equilibrium transitions among these states have 
been observed and studied in heterogenous networks. 

The phase transitions among different collective states 
in heterogenous network exhibit the typical feature of 
multistability, i.e., they may coexist for given parameters 
and depend on the choice of initial conditions. This inter¬ 
esting behavior is closely related to the first-order transi¬ 
tion, especially the recently concerned explosive synchro¬ 
nization. Multistability and discontinuous transitions in¬ 
dicate the competitions of miscellaneous attractors and 
their corresponding basins of attraction in phase space. 
For a network of coupled oscillators, the miscroscopic dis- 
cription of the dynamics of oscillators should be made in 
a high-dimensional phase space, which is very difficult 
to deal with. The key point in understanding macro¬ 
scopic transitions is the projection of the dynamics from 
this high-dimensional space to a much lower-dimensional 
subspace. This can be executed by introducing appropri¬ 
ate order parameters and building their dynamics. Espe¬ 
cially the significant work of m , in that case the origi¬ 
nal infinite-dimensional problem were reduced exactly to 
a flow in OA manifold which was described by several 
coupled ordinary differential equations(ODE). Moreover, 
for the case with finite oscillators, an ensemble method 
based on the OA manifold is also developed. Interest¬ 
ingly, the time-asymptotic attracting behavior of the full 
system and the attractor bifurcations are contained in 
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this reduction and were quantitatively described by the 
ODE formulation. 

Recently the abrupt transition from incoherent state to 
synchronization in network which the frequencies of os¬ 
cillators on nodes are positively correlated to the node’s 
degrees has attracted a lot of attention |5|. This ex¬ 
plosive transition to synchronization has been observed 
numerically on scale-free networks and experimentally in 
electronic circuits [3 [TO,. Moreover, it was also found 
that the first-order transition can be changed and nu¬ 
merous phase transitions can be observed by controlling 
the phase shift |22| . Numerous efforts have been made to 
understand the mechanism of explosive synchronization 
from different viewpoints such as the topological struc¬ 
tures of networks, the coupling functions among nodes, 
and so on jTUHIS]- Although there have been numerous 
discussions and even analytical solution on understand¬ 
ing the explosive synchronization, it is still difficult to get 
an analytical insight in a high-dimensional phase space 
composed of dynamical variables of oscillators, and the 
dynamical analysis of transitions among different syn¬ 
chrony states is still lacking. 

In this paper it is our mission to study the phase 
transitions among oscillators in a star network with the 
Sakaguchi-Kuramoto couplings by considering the ef¬ 
fect of the phase shift among coupled oscillators. The 
star topology is the main topological characteristics of 
many heterogenous networks, e.g., the scale-free net¬ 
works. The dynamics of star networks of oscillators 
is analytically studied by proposing the ensemble or¬ 
der parameter (EOP) and building the equations of mo¬ 
tion of EOP for networks with a finite size,, which 
accomplishes a great reduction from microscopic high- 
dimensional phase dynamics of coupled oscillators to a 
macroscopic low-dimensional EOP dynamics. Based on 
the EOP dynamics, we further reveal, by using tools 
of nonlinear dynamics such as time reversibility anal¬ 
ysis and linear stability analysis, that numerous phase 
transitions among different synchronous states in this 
model. We found three typical processes of the transi- 
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tions to the synchronous state, i.e., the transitions from 
the neutral state, the in phase state or the splay state 
to the synchronous state, and the continuous process 
of de-synchronization together with a group of hybrid 
phase transitions which are discontinuous with no hys¬ 
teresis. Furthermore, we study the synchronization dy¬ 
namics of two coupled star networks and find a two-step 
synchronization transition, which can be macroscopically 
resorted to the dynamics of coupled EOP equations. 

The present paper is organized as follows. In Section 
II, we illustrate the Ott-Antonsen ansatz and the EOP 
equations for network with a finite size. In Section III, we 
introduce the Sakaguchi-Kuramoto model with star net¬ 
work connections and apply the EOP approach to project 
the dynamics of coupled oscillators to the phase space of 
EOP. In Section IV, the properties and mechanics of the 
phase transitions are analyzed in the plane of ensemble of 
order parameter. Moreover, in the last section, the star- 
connected network is also discussed, and a two-step syn¬ 
chronization transition is found and analyzed with simi¬ 
lar methods for the single star network. We believe that 
this is a meaningful step to understand the synchroniza¬ 
tion transition in heterogeneous networks. 

II. OTT-ANTONSEN ANSATZ 

Recently the issue of low-dimensional dynamics was 
reopened in work of jTSj with the discovery of an ansatz 
that collapses the infinite-dimensional Kuramoto model 
to a two-dimensional systems of ODEs. 

To illustrate the ansatz briefly, let us apply it to the 
class of identical oscillators governed by the equation 

^ = fe^ + g + fe~^, j = ( 1 ) 

where / is any smooth, complex-valued 27r-periodic func¬ 
tion of the phases ipi , • • • , tpN and g has to be real valued 
since ipj is real. In the limit N — > oo, the evolution of 
the system (1) is given by the continuity equation 


where r and $ are defined via 

a = re 1 ®. (6) 

In geometrical terms, the ansatz (4) defines a submani¬ 
fold in the infinite-dimensional space of density functions 
p. This Poisson submanifold is two dimensional and is 
parameterized by the complex number a, or equivalently, 
by the polar coordinates r and $. 

The intriguing face discovered in work 118] is that the 
Poisson submanifold is invariant: If the density is initially 
a Poisson kernel, it remains a Poisson kernel for all the 
time. To verify this, we substitute the velocity field (3) 
and the ansatz (4) into the continuity equation (2), and 
find that the amplitude equations for each harmonic e m ^ 
are simultaneously satisfied if and only if a(t) evolves 
according to 

a = i{fa 2 + ga + f). (7) 

This equation can be recast in a more physically mean¬ 
ingful form in terms of the complex order parameter de¬ 
noted by < z> and defined as the centroid of the phases 
t f> regard as points e' 1 ^ on the unit circle: 



By substitute Eq.(4) into Eq.(8) we find that < 2 >= a 
for all states on the Poissan submanifold. Hence, jz l 
satifies the Riccati equation 

< z >= i(f < z > 2 +g < z > +/). (9) 

at 

When / and g are functions of < z > alone, as in mean- 
field models, Eq.(9) constitute a closed two-dimensional 
system for the flow on the Poissan submanifold. 

For the finite K, we can get the similar property of 
the system through an ensemble way. By choosing an 
ensemble consisting of systems with the same parameter, 
we define an ensemble order parameter as 


dp 

dt 



( 2 ) 


where the phase density p(<f>, t) is defined such that 
p{(j),t)d(j) gives the fraction of phases that lie between 
</> and <p + dcj) at time t, and where the velocity field is 
the Eulerian version of Eq.(l), 

u^t) = fe^ + g + fe~^. (3) 

In the work of [ 18 ] , suppose p is of the form 

OO 

p(<M) = 2^{l + 5>(t)"e in * +a(i)"e- i "*)} (4) 


for some unknown function a that is independent of <f>. 
Note that Eq.(4) is just an algebraic rearrangement of 
the usual form for the Poissan kernel 


P{4>) = 


1 


1 - H 


2-7T 1 — 2r cos (</> — <h) + r 2 


( 5 ) 


<«(*)> = (io) 

j=i 

where < • > means ensemble average, K is the number of 
oscillators in a single system and <Pj,j = 1 ,... ,K are the 
phase variables for oscillators. With the approximated 
OA ansatz for network with finite size, if the initial phases 
of oscillators are chosen randomly from a Poisson kernel, 
we have the dynamical equation for the ensemble order 
parameter ( a(t )) as 

= *(/(a(t)> 2 + g(a(t)) + /), (11) 

with assumptions ( fa(t ) 2 ) ~ f(a(t)) 2 and ( ga(t )) ~ 
g(a(t)). With the assumptions work, this equation for 
ensemble order parameter can be used to get the low- 
dimensional collective behaviors of the system with finite 
size. 
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III. STAR NETWORK WITH PHASE SHIFT 


In heterogeneous networks, e.g., the scale-free net¬ 
works, nodes exhibit essentially different linking prop¬ 
erty, among which some nodes with dominantly large de¬ 
gree, i.e., hubs, play a dominate role. The simplest case of 
heterogeneous networks is the star topology, where a hub 
connects with other leaves with unit degree, can grasp 
the essential properties of heterogenity. As our starting 
point of discussion, we take a star-network of coupled 
phase oscillators with phase shift as our working model. 
The effect of phase shift among coupled oscillators has 
been extensively investigated in recent years, while this 
has been seldom discussed in star networks psou. 

Consider coupled phase oscillators with a star connect¬ 
ing topology. The star network is composed of a central 
node (i.e., the hub) and K leaf nodes. Each of the leaf 
nodes connects solely to the hub. Thus, the degree of 
the leaves is hi = 1 (i = 1,... ,I \) while that of the hub 
is kh = K. Suppose that the natural frequency of the 
oscillators are proportional to their degree, the equations 
of motion for the hub and leaf nodes read 

K 

9h=uj h + AV sin (6j - 6 h - a), 

U ( 12 ) 

6j = w + A sin (Oh — 9j — £*)> for 1 < j < K, 

where 9h , 9j and uih, w are instantaneous phases and nat¬ 
ural frequencies of the hub and leaf nodes respectively. A 
is the coupling strength. K is the number of leaf nodes 
connected with this hub and a is the phase shift. It was 
reported that abundant collective dynamics appear in the 
model with phase shift, such as synchrony can decay or 
incoherence can regain stability with increasing coupling 
and multistability between partially synchronized and/or 
the incoherent state can appear in the globally coupled 
network. The star network was original to study the 
characteristics of the explosive synchronization without 
phase shift, to explore the features of the transitions in 
the star network, it’s necessary to investigate the system 
with phase shift. In such conditions, several new states 
are found, the characters and the stable conditions of 
them are analyzed in detail in the following part. 

When K = 1 the model is solvable. With the definition 
of phase difference tp = 9h — 9\ and difference of natural 
frequencies Aw = w/, — w, we have the dynamics for tp as 

tp = Aw — 2A sin(</?) cos(a) (13) 


where cos(a) is only an regulatory factor added to cou¬ 
pling strength A. The synchronous state of the two oscil¬ 
lators means 9h(t) = 9\(t) + c where c is a constant, or 
the equivalent tp = 0. Hence the existence condition for 
synchronization reads 


A > A c 


Aw 

2cos(a) 


(14) 


unless cos(a) = 0. With A > A c , the constant c can 
be solved as tpo which satisfies the condition tp = 0, as 


tpl = arcsin(Aw/(2Acos(a))) or tp^ = 7r — tp J. The two 
possible synchronous states is either stable or unstable 
which determined by — cos(a) cos(</?o) by the linear sta¬ 
bility analysis. When cos(a) > 0, the state tp J is stable 
and tp% unstable. When cos(a) < 0, the case is just op¬ 
posite. In the case of cos(a) = 0, we have tp = Aw which 
is a constant and the system cannot be synchronous un¬ 
der any circumstances. Obviously, the phase shift a has 
great influence on the dynamics of the model in this case. 

For the cases K > 2, with phase differences Pj = 
Oh — 9j, the star network with phase dynamics can be 
transformed to an all-connected network with phase dif¬ 
ference dynamics as 

K 

Pi = Aw — A sin(<^ + a) — A sin(<p. i — a), 1 < j < K 
i=i 

(15) 

To describe the synchronous state or non-synchronous 
state, we define the order parameter as 

1 K 

z(t) = r(ty m = — ^ ei{Vj) - ( 16 ) 

A 3=1 

It is instructive to rewrite Eq.(13) as 

<Pi = +g + fe -**, j = 1, • • ■ , K, (17) 

where i denotes the imaginary unit and / = i— , g = 
Aw — A Kr sin($ + a). 

For the finite K, we can get property of the system 
through an ensemble way. By choosing an ensemble con¬ 
sisting of systems with the same parameter, we define an 
ensemble order parameter as 

z{t) = <r(t)e 4 *W> = <1 e iM ) = (Z(t )). (18) 

3=1 

With the approximated OA ansatz for network 

with finite size, if the initial phases of oscillators are cho¬ 
sen randomly in an fixed interval, we have the dynamical 
equation for the ensemble average of order parameter z(t) 
as 

z= -^e~ ia z 2 +i(Auj-\Krsm{$ + a))z + ^e ia , (19) 

and Z(t ) ss z(t) from which we can use z(t) to determine 
the state of the system instead of Z(t). Eq.(8) describes 
the collective dynamics of Eq.(l) in terms of the order 
parameter, different solutions of Eq.(8) build correspon¬ 
dences with diverse collective states. Following we will 
analyze the characters and stable regions of those collec¬ 
tive states with different parameters. 


A. Steady solutions 

The ensemble order parameter equation describes the 
collective dynamics of the star network. According to the 
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Eq.(8), two steady solutions are found. The characters 
and the stable region in the phase diagram of them are 
analyzed as follows. 

With Eq. (19), setting z = x + iy, we can describe the 
system in a x-y plane as 


x = A(- + K) cosay 2 — — cos ax 2 


+ A (K — 1) sin axy — Aw y + — cos a, 
y = A(- — K) sin a a; 2 — — sin ay 2 
— A (K + 1) cos axy + Aw x + — sin a, 


( 20 ) 


which can be studied analytically. 

By setting x = 0 and y = 0, we have the four fixed 
points noted by ( Xi,yt ) with 


— sin a Aw ± A sin a 
\{2K cos 2a + 1) 

sina S -^B ±K{ S -^B-sm2a 2 ) 

A Asina(/\ 2 + 2cos(2a)/v + 1) 

— cos a Aw ± A cos a 
\{2K cos 2a + 1) 

—Aw(— cos a ± sin aB — K cos a) 

A(Jv 2 + 2 cos(2a)A' + 1) 


where A = \J—2K\ 2 cos 2a — A 2 + Aw 2 and B = 
y/ A 2 T A' 2 A 2 + 2KX 2 cos 2a — Aw 2 which also give the 
existent condition for the fixed points as for point 1 and 
2, as 


A < Ai 


Aw 

\J2K cos 2a + 1 ’ 


( 22 ) 


and for point 3 and 4, as 


A > A 2 


Aw 

\JK 2 + 2K cos 2a + 1 


(23) 


For point 3 and 4, there is an additional relation x 2 +y 2 = 
1 while for point 1 and 2, x 2 + y 2 may be greater or lower 
than 1. For the definition of z, the point is allowable to 
describe the state of the system if and only if x 2 +y 2 < 1. 

As the analysis we did for the simple case with K = 
1, the linear stability analysis can be used for the fixed 
points, with which the Jacobian matrix is calculated, 

T _ / J 11 Jl2 \ 

~ \Jll J 21 ) 


with Jn = — A cos air + A (K — 1) sin ay, J 12 = A(1 + 
2K) cos ay + A (K — 1) sin ax — Aw, J 21 = A(1 — 
2K)s'max — A (K + 1) cos ay + Aw, J 22 = —A sin ay — 
A {K + 1) cosax ,the eigenvalues are 


/? 1,2 


J\\ + J22 i + J22) 2 — 4(Jn J22 — J12J21) 

2 ' 
(24) 



FIG. 1: (color online) (a) The time evolution of shnp^t) with 
a = — 0.4-7T, A = 2, ( b ) The order parameter against the cou¬ 
pling strength with different initial states for a = 0.0, (c) The 
time evolution of sin0i(t) with a = 0.l7r,A = 0.5, ( d ) The 
lyapunov exponents of the network with a = 0.17T. The size 
of the star network is N = 11. 


Then the stability of the fixed points can be judged from 
the eigenvalues of the matrix. If and only if the eigenval¬ 
ues both are negative ,the fixed point is stable. 

For point l,the stability condition reads as 


\<4 


Au 

y/2K cos 2a + 1 


for a £ (a Q , 0) 


(25) 


where Oq = — arccos(— l/K)/2 and when —n/2 < a < 
a(j" this point is always stable. It is worthy to note that 
the stability condition for it is the same with existence 
condition in Eq. (22) for a < 0. As for a > 0 it is always 


unstable. In the critical case with a = 0, the point has 
the neutral stability. 

For point 2, the stability condition reads as 

-Jora £ («+, tt/2), (26) 

K cos 2a + 1 

where a(j" = arccos(— l/K)/2. 

For point 3,the stability condition reads as 



Aw 


K cos 2a + 1 ’ 
—Aw 


K cos 2a + 1 ’ 


fora e (a 0 ,0), 
fora G (a^,7r/2), 


(27) 

and when 0 < a < aj the point is always stable. It 
is worthy to note that, the critical coupling A)j" c is spe¬ 
cial that when A > A+ the point 2 is stable while point 
3 is unstable, and when A < A+ the case is just oppo¬ 
site. At last, for point 4, it is always unstable under any 
circumstances. 

The fixed points of order parameter correspond to the 
stable solutions of coupled oscillators. For points 3 and 
4, with 1 2 1 = -y/x 2 + y 2 = 1 these two points correspond 
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FIG. 2: (color online) Phase diagram of system. Regime SS 
is stable synchronization region. Regimes SPS and IPS are 
stable region for splay state and in phase state respectively. 
The stable region for neutral state is too narrow to plot with 
only a = 0, if. The coexistence regime of splay state and 
synchronous state is plotted by shadow. 


to the synchronous state (SS) of the system which means 
all the phase differences between hub and leaf nodes are 
the same and fixed as = constant for 1 < j < K. 

Hence the point 4 is always the unstable synchronous 
state and point 3 is the stable synchronous state in some 
special regions as we showed. 

For points 1 and 2, if \z\ = x 1 + y 2 > 1 then these 
states is not attainable because the order parameter a 
is bounded by \z\ < 1. If |z| < 1 then the points are 
related to the splay state (SPS) which means all the phase 
differences satisfy a single function as <Pj(t) = ip(t + ^) 
with T the period of ip(t) for 1 < j < K , as shown in Fig. 
HH For — 1/27T < a < 0, the point 1 is the stable splay 
state and point 2 is the unstable one in same regions and 
for 0 < a < 1/27T the case is opposite, the stable splay 
state corresponds to the point 2. 

All the existent and stable region for the synchronous 
state and splay state consist the phase diagram plotted 
in Fig. 2. It can be described by the four regions of the 
phase shift a. When —l/2n < a < a q , the splay state 
corresponding to point 1 always exist and is stable and 
the synchronous state exist if A > X ec = X 2 but always 
unstable. When < a < 0, the splay state exist and 
is stable with A < X{ = Ai, and the synchronous state 
exist with A > X ec but unstable unless A > Xj c > X ec . 
It is noted that there is a coexist region for the splay 
state and synchronous state when Aj c < X < X[. In 
the third region, when 0 < a < a,}, the splay state is 
always unstable which only exist for A < Ai, and the 
synchronous state exist and is stable for A > A ec . At 
last, when < a < 1/27T, the splay state exist all the 
time but only stable when A > A+, and the synchronous 
state only exist and stable in the region A ec < A < A+, in 
another hand. 



FIG. 3: (color online) (a) The forward and backward contin¬ 
uation diagrams in different initial conditions with a = 0, (6) 
The forward continuation diagram with a = ±0.l7r. 


B. Periodical solutions 

With the analysis of fixed points, there are still two 
regions in the phase diagram where the stable state of 
system unclear, the one with 0 < a < 1/27T and A < A ec 
where only exist unstable splay state, and the ones with 
a = 0,±1/27t which is the critical cases. It is clearly 
that the periodical solutions must exist in those regions. 
Following, we will find the periodical solution and analyze 
the character and stable region of them. 

For the first case, when A < A ec and 0 < a < 7t/ 2, it 
is shown that there is no stable fixed point exist. How¬ 
ever, from the numerical simulation and lyapunov expo¬ 
nents for the trajectories, shown in Fig. 1(d), we find the 
largest lyapunov exponent is zero and the others are neg¬ 
ative which means there exist a limit cycle as the stable 
state and the volume of phase space is convergent. 

To find the limit cycle, we transfer the Eq .(8) from 
x-y plane to polar coordinates by setting z = re l4> , 

r = — ^(r 2 — 1) cos(<f> + a), 

$ = —^(r H—) sin('I> — a) + Aw — XKr sin($ + a). 

2 r 

(28) 

It is easy to see that there is a limit cycle solution with 
r = 1 and periodic phase <h(t) for A < A ec . This state we 
call in phase state(IPS), where all the phase differences 
are the same as tpj = ip(t ) for 1 < j < K. This state 
means in the dynamical model Eq.(l) all the phases of 
leaves nodes are identical, as shown in Fig. [ljc) . In this 
case, the model Eq.(l) is reduced to the case with K = 1 
and the stability of the state can be got by the Flouquet 
theory for limit cycle as it is stable for 0 < a < 7t/ 2 and 
unstable for — 1 / 27 T < a < 0. 

For the second case, when a = 0, ±l/27r, these are the 
critical cases for the system. In the above analysis for 
K = 1, with a = ±1/27r we have the conclusion that 
the system will never be synchronous. Meanwhile, the 
fixed point 1 is found to be neutral stable when a = 
0, ±1/27t. With numerical simulation, we find that there 
is a large class of states in the critical cases with the 
order parameter can be almost every number between 0 
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FIG. 4: (color online) Phase plane of Eqs.(29l with A u = 9, 
K = 10, a = 0 (a) with A = 1.5, ( b ) with A = 1.9. Red lines 
are x = 0, and green lines are y = 0. The intersections of 
x = 0 and y = 0 are fixed points A, B, C, D. Trajectories 
with different initial values are marked by 


and 1 , determined by the initial situation of the nodes, 
shown in Fig. 1(b). When a = 0 there are transitions 
from the special state to synchronous state at different 
coupling strength. And when a = ±l/27r the system 
will always in the states whatever the coupling strength 
is. Moreover, the sum of Lyapunov exponents for these 
states is zero which means the volume of phase space 
is conserved and hence we call them the neutral states 
(NS) . 

To understand the mechanics of these neutral states, 
in the framework of approximated Ott-Andeson ansatz 
and the plane of ensemble order parameter, we get the 
dynamical equation in with a = 0 as 


is called reversible trajectory. For any reversible trajec¬ 
tory, the Lyapunov exponents have the sign-symmetry 
form and the volume of phase space in the vanity of it 
are conserved as we discovered with numerical simula¬ 
tion. 

For our planer system bounded by unit circle with the 
invariant set as x = 0 , y > 0 , when the only fixed point 
is neutral stable, every trajectory is closed periodic one 
which belongs to the neutral states. This is what happens 
when a = ± 1 / 27 t for all the coupling strength and a = 0 
with the region A < X ec . For a = 0 and A > A ec , there is 
an coexist region for the synchronous state and neutral 
states as the critical cases for the coexist region for splay 
state and synchronous state. 

Up to now, we get the phase diagram completely as 
shown in Fig .2. There are three main states, namely, 
the splay state, in phase state and synchronous state 
and one kind of critical state, the neutral state. The 
stable region for each state is plotted in the phase dia¬ 
gram with the boundaries we get analytically above both 
from existence and stability conditions. There is also a 
coexistence region for splay state and synchronous state, 
with the critical case for neutral state and synchronous 
state. The variety of states in the phase diagram leads 
to various phase transition of the system. 


C. Processes of synchronization and 
de-synchronization 


x = X (K +-)y 2 — —x 2 — A uy + —, 

v 2 y 2 y 2’ (29) 

y = —X(K + 1 )xy + Acox. 

By setting x = 0 and y = 0,when A < X ec only one fixed 
point exist inside the unit cycle in the plane as the point 
1 we analyzed. 

According to the linear stability the fixed point is neu¬ 
tral stable. It is worthy to note that if we define an time 
reversal transformation as R : (t,x,y) (—t,—x,y), 

with R effect on both sides of the equations, the dy¬ 
namical equations remains invariant and hence they are 
called the time reversal dynamical system. This is a spe¬ 
cial kind of symmetry of dynamical system, which is also 
called the quasi-Hamiltonian m- This symmetry gives 
the system many interesting properties. 

Note that, the time reversal transformation R can be 
resolved into R = TW with T : t >->■ —t and W : (x, y) e-i 
(—x, y). Hence the invariant set for W is the y axis with 
x = 0,y > 0. For any trajectory cross this invariant set, 
according to the time reversal symmetry, the forward tra¬ 
jectory and the backward trajectory are symmetrical. If 
the forward trajectory evolves to an attractor then there 
will be an symmetrical repeller of the system, which is 
the same attractor for backward trajectory. If the tra¬ 
jectory cross the invariant set more than once, then the 
forward and backward trajectory will coincide with each 
other, forming the periodic solution for the system, which 


As we see in the phase diagram, the model we 
use has great variety of phase transitions among the 
four special states, as the transitions to synchronous 
state, namely synchronization, and correspondingly de- 
synchronization, also the transition between the same 
kind of states when crossing the critical state with a = 
0, ± 57 r. Some of the states have coexistent regions, which 
lead to abrupt transitions among them and hysteresis be¬ 
haviors, while the others lead to transitions either abrupt 
or continuous, according to the kind of bifurcation at the 
critical point. 

Among all the transitions, the processes of synchro¬ 
nization are most important. In the phase diagram it is 
clear that there are three processes to synchronous state 
from splay state, in phase state and neutral state respec¬ 
tively, with fixed phase shift a and increasing coupling 
strength A as usually assumed. 

The first, we take an view on the synchronization pro¬ 
cess as neutral state to synchronous state with a = 0 . 
The synchronization process is discontinuous known as 
the explosive synchronization which has attracted much 
attention recently. According to numerical computations 
this kind of transition is revealed abrupt, and there is a 
hysteretic behavior at the onset of synchronization. X b c 
and X{. are the backward and forward critical coupling 
strengths respectively, where X b = X 2 and X{ depends on 
initial states as shown in Fig. [3|a). The upper limit of 
X{ is denoted by X[. As A > X(, the synchronization 
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FIG. 5: (color online) Phase plane for Ao; = 9, K = 10, 
a = —0.17r, (a) with A = 1.8, (b) with A = 2.17. Red lines 
are x = 0, and green lines are y = 0. The intersections of 
x = 0 and y = 0 are fixed points A,B,C,D. Trajectories with 
different initial values are marked as V. 


FIG. 6: (color online) Phase plane for A uj = 9, K = 10, 
a = 0.3-7T, (a) with A = 0.5, ( b ) with A = 1.5. Red lines 
are x = 0, and green lines are y = 0. The intersections of 
x = 0 and y = 0 are fixed points A,B,C,D. Trajectories with 
different initial values are marked as V. 


state is globally attractive. Apparently it is difficult to 
understand the process of it according to the self consis¬ 
tent method, especially for the hysteresis behavior and 
coexistent region. 

However, as we see, through the method of OA anstaz, 
the coexistence of neutral state and synchronous state is 
pronounced, and the basins of attraction of each state 
can be determined also. For the time reversal symmetry 
the system has, we have defined an time reversal trans¬ 
formation as R = TW : ( t , x , y) (— t, —x, y) and the 
invariant set for W is x = 0, y > 0. In the case with the 
planer system bounded by unit circle, every trajectory 
either evolves to an stable point in the limit cycle or is 
closed periodic one with neutral stability. The separation 
between the two states is the condition for the trajectory 
to evolve only one round, as 


dx 

dt 


d 2 x 

= 0, dt? X> ° 


dy_ 

dt 


=°^y >0 


(30) 


The trajectories are plotted in Fig. |4j The separation 
is the up red line with which the basin of attraction for 
neutral state and stable synchronous state are separated. 
The points D and B are the points 1 and 2 corresponding 
to the splay state and points A and C are the points 3 and 
4 corresponding to the synchronous state. When a = 0, 
point D is neutral stable and point B is an saddle. Points 
A and C are linear stable and unstable respectively. The 
neutral states are the closed periodic orbits around the 
point D. 

With the changing of parameters as coupling strength 
A or difference of natural frequency Aw the separation 
line changes and phase transition will happens between 
these two steady states with a single initial state with 
marks * in Fig. [4] The transition is abrupt or explosive 
apparently, as shown in Fig. [3ja) , with hysteresis. And it 
is easy to take the basion for neutral state approximated 
by the circle which has its center in point D and radius 
as the length of line B — D. 

When a < 0 it can be shown that the point D becomes 
linear stable. If it is the only stable point in the phase 


space, all the trajectories will evolve to it as splay state of 
the system. In the coexistent region with another stable 
point A, which corresponding to the synchronous state, 
both the two stable points has their own basion of at¬ 
traction. In this case, the point B plays the role of the 
separations of the attraction regions for both stable state. 
It is shown in Fig. [5] With the parameter changing, the 
point B as the separation changes. The transition from 
splay state to synchronous state is also an abrupt one, 
shown in Fig. §b), also with hysteresis. 

When a > 0 the points D and B are both unstable and 
before the appears of points A and C there is no stable 
points in the unit circle. The trajectories in this case will 
all evolves to the boundary which is a limit circle for the 
planar system. Hence when a stable point A and unstable 
point C appears in the unit circle the transition from the 
in phase state to synchronous state is continuous as the 
ghost effect of stable fixed points which is shown in Fig. 
[3jb) . The process of phase transition in the planer phase 
space is shown in Fig. [6j 

The three ways to synchronization discussed above are 
either abrupt or continuous. As we see, the different 
routes to synchronization is based on the different kinds 
of incoherent state, namely, the neutral state with a = 0, 
the splay state with a < 0 and the in phase state with 
a > 0. The small change of a around zero can change 
the process of synchronization dramatically. This is the 
reason we are interested in the model with phase shift. 

Contrary to the conventional belief that the system 
will always be synchronous if we set the coupling strength 
large enough, the synchronous state is unstable when A > 
A+. It is clear in the phase diagram. The process is de- 
synchronization which is shown in Fig. [7](a) and (b). The 
order parameter decrease from 1 and effective frequencies 
of hub and leaf nodes are divided at the same time. This 
a transition from synchronous state to splay state. 

When A < A+, the system is in synchronous state, 
namely the point A in the phase plane. At the critical 
case A = A+, the red line of x = 0 and the black line of 
y = 0 are tangent in the unit cycle in fixed point A where 
the other fixed points D and C are still unstable. If we 
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FIG. 7: (color online) (a) The order parameter against the 
coupling strength with a = 0.3-7T. (b) The effective frequency 
of the hub nodes(the green) and the leaf nodes(the yellow) 
against the coupling strength with a = 0.37T. Phase plane 
for Aw = 9, K = 10, a = 0.37T, (c) with A = 3, (d) with 
A = 5. Red lines are x = 0, and green lines are y = 0. The 
intersections of x = 0 and y = 0 are fixed points A,B,C,D. 
Trajectories with different initial values are marked as V. 

increase the coupling strength and let it larger than A+, 
the two nullclines will intersect in two fixed points A and 
another one, the point B. The point A loses its stability 
and the new point B is a stable one corresponding to the 
splay state. The process from the synchronous state to 
splay state is finished by this bifurcation continuously. If 
we reverse this process by decreasing coupling strength, 
it is not doubt a process of synchronization, although we 
usually not see it in this view. 


D. Hybrid transition 

Up to now, we have discussed four kinds of phase tran¬ 
sitions, either abrupt with hysteresis or continuous. All 
of them are concerned with synchronous state and we 
follow the tradition with fixed phase shift a and increas¬ 
ing coupling strength A. However, there is still a class 
of phase transitions which seems to be interesting for us, 
with fixed coupling strength A and changing the phase 
shift a, which allows us to think the process across the 
critical lines with a = 0, ±1/27t. 

Surprisingly, at a glaze, this crossing process always 
leads to a kind of hybrid transition which is abrupt but 
without hysteresis, from splay state to in phase state, 
and even between synchronous states and splay states in 
the different sides of these lines, as shown in Fig|8ja) for 
transitions between splay states. 

However the reasons of it is not very difficult to un¬ 
derstand. Based on the analysis above for the cases with 
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FIG. 8: (color online) (a) The order parameter against time 
in different phase shift a with A = 5. Phase plane for A uj = 9, 
K = 10, A = 5, ( b ) with a = 0.49-7T, (c) with a = 0.50-7T, (d) 
with q = 0.517T. Red lines are x = 0, and green lines are 
y = 0. The intersections of x = 0 and y = 0 are fixed points 
A,B,C,D. Trajectories with different initial values are marked 
as V. 


a = 0, ±l/27r, they are protected by the time-reversal 
symmetry as the critical cases. Once the phase shift 
crossing the critical cases, all the stability of the solu¬ 
tions will turn to the opposite, which is typical for hybrid 
transition. 

Take the phase transition between splay states as an 
example. In the phase plane of order parameter, there 
exist four fixed points for each side of the system. One 
stable point and one unstable point inside the unite cy¬ 
cle with two saddle points in it. As shown in Fig{8](b), 
the point D is stable and B unstable. When crossing the 
critical case, both the points D and B are neutral stable 
in Figgc). And in other side of the critical case, the 
stability of them turn to the opposite which means the 
saddle points are still the saddle points but the stable 
point becomes unstable and vice versa. It is shown in 
Figj8|d) with the point B the stable one and D the un¬ 
stable one. The process from (b)-(c)-(d) or (d)-(c)-(b) in 
Figj8] is a typical hybrid transition. 

Noting that all the states are the collective behaviors 
of the phase oscillators in the star network, the phase 
transitions here are not the trivial ones at all. This is 
largely related to the symmetry of the system, both the 
identity of the oscillators which makes the approximate 
OA ansatz appropriate and the time-reversal symmetry 
for which phase shift affects mostly. Among all the cases, 
the most important one is the case with a = 0, which has 
the time-reversal symmetry and also the original model 
to be used for discussion of synchronization. 
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IV. STAR-CONNECTED NETWORK 

In the heterogeneous network especially scale-free net¬ 
work, hubs play the dominate roles. The single star net¬ 
work which get the essential properties of hubs is very 
useful to explore the main behavior of heterogeneous 
network which we have discussed above. However, the 
heterogeneous network such as scale-free network have 
several hubs, and hubs has influences on each other. To 
better describe the phase transitions in the heterogeneous 
network, we can combine several star networks together 
and let the hubs connected directly, namely the star- 
connected network. 


A. Model and the simplification 

Considering the simplest network composed by two 
star networks, namely star 1 and star 2. The model reads 
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K i 

0/n = w/n +Ai ^ sin (6^ -9 hl ) + X 0 sm(e h2 - 9 hl ), 

■ 7*1 = 1 

6 h = u/r + Ai sin (6 hl - 9 h ), for 1 < ji < Ad, 
k 2 

0h 2 = Uh 2 + A 2 ^2 sin (9 h -9 h2 ) + A 0 sin (9 hl - 0/ l2 ), 

32 = 1 

9j 2 = w 2 + A 2 sin (i 9 h2 - 0 i2 ), for 1 < j 2 < K 2 , 

(31) 

where 9h 1 ,9h 2 ,9j 1 ,9j 2 and Uf ll ,Uh 2 iUj 1 ,Uj 2 are phases 
and natural frequencies of the hub and leaf nodes in star 
1 and star 2 respectively. K\,K 2 and Ai, A 2 are the num¬ 
bers of leaf nodes and coupling strength for star 1 and 2. 
Ao is the coupling strength between the two stars. 

With a = 0, in the single star network without phase 
shift, the phase transition is from neutral state to syn¬ 
chronous state with abrupt transition. However, when 
coupled together, setting Ai = A 2 = Ao = A, the syn¬ 
chronization process with increasing A for a single star 
network is divided into two parts. A intermediate state is 
observed that the leaf nodes of each star will synchronize 
before the synchronous state appears. The process ap¬ 
pears to be from neutral state to in phase state abruptly 
then to synchronous state continuously, if we borrow 
the definition for the single star network, as shown in 
Fig^a). 

The second step of transition from in phase state to 
synchronous state is also observed for the case with a 
single star. This can be seen as a bifurcation as the 
synchronous state appearing in the in phase state and 
is a consequence of the inside coupling strength Ai and 
A 2 large enough. On another hand, if we put the inside 
coupling strength Ai = A 2 = A fixed, and changing the 
coupling strength Ao between the stars, there is only the 
first part of synchronization, the abrupt transition to the 
state with leaf nodes synchronized, as shown in Fig[9](b). 
Hence the two parts of synchronization can be seen as 


FIG. 9: (color online) The order parameter of star 1 n and the 
local order parameter of leaf nodes of star 1 tli against the 
coupling strength A or Ao with different coupled forms when 
K i = 10, K 2 = 9, Ai = A 2 = A. (a) Ao = A, bidirectional 
coupled with the star 2. ( b) A = 1.3, bidirectional coupled 
with the star 2. (c) A = 1.3, unidirectional coupled by the 
star 2, (d) A = 1.3, unidirectional coupled by the hub of star 
2 . 


two kinds of mechanics, caused by Ai and A 2 , or Ao, and 
it is possible to study the two steps of phase transition 
through controlling Ay 2 or Ao respectively. 

For each star, the influence of the other is only a term 
added to the dynamics of the hub, which is small com¬ 
pared to the total influence of its leaf nodes. Hence, 
before the transition to the state with leaf nodes syn¬ 
chronized, the state of each of the stars is almost the 
same with neutral state for a single star network, which 
is also can be seen in Fig[9](b) before the transition. To 
find the mechanics of the first transition from neutral 
state to in phase state, we can assume the star 2 is just 
not affected by star 1 at all, consequently in its neutral 
state, and focus on how the state of star 1 is affected 
by it through the unidirectional coupling. As a result, 
with the increasing of Ao, the appearance of the phase 
transition in star 1 is still observed as shown in Figj9|c). 


Furthermore, a more simplified model can be consid¬ 
ered, as the limit of neutral state, with A 2 = 0 that the 
star 1 only driven by the hub of star 2. This simplified 
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FIG. 10: (color online) (a) and (c) The phase space with 
A = 1.3, Ao = 0 or 0.64 respectively. The blue lines are the 
trajectories of the simplification model with the transform 
X(t ) = (3 + x{t) * cos(a(t))), Y(t) = (3 + x(t ) * sin(a(t))), 
Z(t) = a(t), other color lines are the boundary of it. (6) and 
(d) The Poincare section of them when a = ■ The red 

points are the trajectories in this section, the blue lines are 
x = 0. 


model reads 


K i 

d hl = u} hl + X ^2 sin ( 9 h - 9 hl ) + A 0 sin(^ 2 - 9 hl ), 

jl = l 


9 h = wi + A sin (0 hl - 0 h ), for 1 < ji < K 1 , 


as 


zi = iAujz i - ^(~i 2 - 1) - ^TjrCzi - zi*)zi - iXoZi sin a, 

• a / NX r \ • 

a = ialu --(zi — z i jzi — Ao sm a, 

2 1 

(33) 

where Aw = Wi — = Uh 2 — uiht ■ 

By setting z\ = x + iy, the system can be described as 


x = -Atxy - ^(x 2 - y 2 - 1) + ^ y 2 + A 0 sinm/, 

?/ = Awa: — Aa;y — NXxy — Ao sin ax , 
d = Au/ — A r Ay — Ao sin a, 

( 34 ) 

It is worthy to note that this new model is quasi- 
Hamiltonian also, with the time reversal transformation 
R which can be resolved into R = TW with T : t \—> —t 
and W : (x,y,a) i —> (—x,y,±n — a). Then the tra¬ 
jectories of the system are also time reversible and the 
volume of phase plane is conserved for reversible trajec¬ 
tories which cross the invariant set a; = 0,a = ±^7r more 
than once. But the system here is three-dimensional, 
which is more complicated than the two-dimensional one 
we discussed above. For instance, around the reversible 
trajectories there can be chaotic trajectories which also 
has conserved volume, similar as the case of Hamilton 
system, and the attractors in this system will not be lim¬ 
ited by fixed point and limit cycle that limit torus and 
chaotic attractor are also possible. To fully analyze the 
system is out of the scope of this paper, we will only fo¬ 
cus on the major property of it related to the transition 
process. 


B. States and transitions 


9h 2 ^h 2 ■ 

(32) 

The drive effect between the star 1 and the hub of star 
2 is depend on the coupling strength Ao- To check the 
rational of the new simple model we still set the A among 
the star 1 unchanged, and reinforce the coupling strength 
Ao from 0. Amazingly to some extent, the abrupt phase 
transition to in phase state in star 1 occurs if Ao achieve a 
critical strength as shown in Figj9](d), just as the original 
model. 

Through the analysis above, it is reasonable to us that 
the major mechanism that responsible for the first part of 
phase transition is maintained in this process of simplifi¬ 
cation. Fortunately, the simplified system is now within 
our ability to analyze analytically through some method. 

To study the mechanism of the phase transition from 
neutral state to in phase state, we define the order pa¬ 
rameter of star 1 as Z\ = rie 1 ® = where 

the phase differences 4>i = — 9 ^, a = 9h 2 — 9h 1 and 

N = K \, with the OA anasatz we can rewrite the model 


By setting x = 0,y = 0 and a = 0 when Ao ^ 0, we 
have several fixed points in the phase space. However, 
they are both nonphysical due to the values of x 2 + y 2 of 
them are greater than 1 which obey the inner constraints 
of the order parameter z that the point is allowable to de¬ 
scribe the state of system if and only if x 2 + y 2 < 1. This 
means that the system can not reach the synchronous 
state as shown in Figj9])d) . 

In another hand, for this system, the state with all the 
leaf nodes synchronized is no doubt a limit circle too, for 
only one degree of freedom left, shown in FigfTl|a). As 
for the neutral stable state, when Ao = 0 it is either 
the periodic trajectory or quasi-periodic torus, shown 
in FigfTo^a), and when A 0 > 0 it is disturbed by the 
coupling between stars which may lead the state to be 
chaotic as shown in in Kig[To|'c). Hence the transition 
from the neutral state to in phase state in this system 
is described by the transition from the chaotic or quasi- 
perotic torus to a limit circle. 

Even though the system has many dynamical struc¬ 
tures as chaotic attractors and so on, what we focus on 
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is only the transition from these two kinds of state which 
is the major factor to the problem. As to the case with 
two dimension dynamics, the abrupt transition and hys¬ 
teresis can be described by the separation line for the 
two states in the coexistence region. And to a system 
with time-reversal symmetry the separation line is right 
the nullclines which prevent the trajectory crossing the 
invariant set twice, as x = 0 for a single star. For the 
three-dimensional system here, the cases are more com¬ 
plicated, but the properties related to the time-reversal 
symmetry is quite similar. The separation for the two 
states, volume conserved trajectory and limit circle, are 
the surface with x = 0. When Ao = 0, it reduce to the 
case for a single star. The surface reads 

— A toy — x 2 — y 2 — 1) + NXy 2 + Aosinay = 0 (35) 

the surface is wavy and it will be the peak or bottom 
when a = f or 4^, which means that the largest distance 
between a point in the surface and a point in the tours 
around center is with a = and the smallest with 4^. 
It is coincide with symmetry analysis that the invariant 
sets are exactly with a = | or 4y. 

For case for a single star, the neutral state depends 
on the initial state and if its trajectory hit the separa¬ 
tion line it will evolve to synchronous state, otherwise it 
will stay as the neutral state. The position of separa¬ 
tion line depends on parameters which gives the basion 
of attraction of synchronous sate and also the critical 
coupling strength for fixed initial condition. For the case 
with three dimension, it is quite similar, except that the 
synchronous state for a single star is now the limit circle 
and the separation line is the separation surface for the 
system. 

For the wavy structure of the separation surface, we 
can only consider the bottom where the surface is closest 
to the torus. For 4f, the variable of x 2 can be formed by 
y that 

x 2 =y 2 + l + j(-\ 0 -Au>)y + 2Ny 2 (36) 

then the range of y is [—1, b with b = 

j (—Ao — Aw). Together with range of y and the rela¬ 
tion with x and y we could get the distance between the 
separation surface and the torus with order parameter ro 


d =r 


min 


- r o 


2 [b 2 + byjb 2 - 4(2TV + 1) - (4AT + 2) {AN + 4)] 
(4 N + 2) 2 


(37) 

If d < 0 the separation surface and the torus have inter¬ 
section and if d > 0 the two surfaces are separated. The 
critical coupling strength for a fixed initial condition is 
the case with d = 0. If the number of N is large enough, 
fmin approximate equal to gx?, which means that when 
the phase transition occurs, the order parameter r has 
a linear relation with the critical coupling strength Aq c . 


r o 



FIG. 11: (color online) (a) The phase space of the simplifi¬ 
cation model with A = 1.3, Ao = 0.9. The blue lines are the 
trajectories with the transform X(t) = (3 + x(t ) * cos(a(t))), 
Y(t) = (3 + x(t) * sin(a(t))), Z{t) = a(t), other color lines 
are the boundary of it. (6) The Poincare section with a = 
of the phase space, the red points are the trajectories in this 
section, the blue lines are x = 0, the fixed point A is the 
intersection of down blue line x = 0 and the unit cycle, (c)- 
( d ) The local order parameter of star 1 against the time with 
different Ao- 


Increasing Ao, r m j n decease and more torus will intersect 
with separation surface and more trajectories will evolve 
to the limit circle. 

As an example and illustration, the process is shown 
in FigjTo] and FigjTT] When A 0 = 0, the star 2 has no 
influence on the star 1, as in Figjl0|a) and (b), and the 
star 1 is in neutral stable state which we have discussed 
above with the trajectories in the phase space regular 
and periodical. When A 0 = 0.7, the star 2 has greater 
influence on the star 1, the trajectories in the phase space 
become chaotic and irregular, from the view of Poincare 
section, the red points jitter and expand greatly. Mean¬ 
while, the separation line, as the down blue one in the 
figure, upwards all the time with the increasing of Ao, and 
the volume of the basin of attraction of in phase state is 
increasing. When the coupling strength Ao between hubs 
is big enough, the trajectories will hit the separation line 
and evolve to the limit cycle which corresponds to the 
in phase state, as shown in FigjTT] when A 0 = 0.9. The 
phase transition occurs now, and all the leaf nodes syn¬ 
chronize in this time. 

Compared with the case with a single star, this dynam¬ 
ical process is quite similar, except three points following. 
First, the fixed point as the synchronous state is related 
to the limit cycle solution as the in phase state. Second, 
the neutral state is now chaotic and inflates because of 
chaos. And third, the coupling between the stars will 
enlarge the basion of attraction of the in phase state, re¬ 
ducing the region for neutral state. Because of the three 
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FIG. 12: (color online) The order parameter of star 1 against 
the coupling strength Ao with different initial states (a) with 
Ki = 10, A '2 = 10, Ai = A 2 = A = 0.5, unidirectional coupled 
by the star 2, (6) with K\ = 10, K 2 = 8, Ai = A 2 = A = 0.8, 
bidirectional coupled with the star 2. 


factors, we have the mechanism of the first step of tran¬ 
sition from neutral state to in phase state, which can 
be checked by the linear relation between Ao c and initial 
condition ?'o, shown in Fig |12| for both the unidirectional 
and bidirectional coupled star networks. 

Regardless of the coupling strength between the two 
star networks either unidirectional or bidirectional, the 
linear relation between the initial order parameter and 
the critical coupling strength Ao c are both observed. The 
linear relations in those models are in agreement with 
the result we get, supporting the analysis above and the 
mechanism we get from the simplified model. 


erogenous network. According to the approximated OA 
ansatz method, we propose the ensemble order parameter 
equation. Based on this equation, the high-dimensional 
phase space can be reduced to a low-dimensional order 
parameter space without additional approximation, then 
the essential dynamical mechanism can be analyzed an¬ 
alytically. In the star network, four characteristic states 
are found by the ensemble order parameter equation and 
the characters and stable regions of them in the phase 
diagram are determined by the time reversibility anal¬ 
ysis and linear stability analysis. Under the action of 
phase shift and coupling strength, extensive phase tran¬ 
sitions among those states are found. Three processes 
of synchronization from splay state, in phase state and 
neutral state to synchronous state respectively, the de- 
synchronization from synchronous state to splay state 
and a group of hybrid transitions in the synchronous 
state, the splay state and from in phase state to splay 
state or the opposite. The processes and mechanisms 
of them are discussed by means of the plane of ensem¬ 
ble order parameter, with planar dynamical of it, all the 
properties of phase transitions can be determined ana¬ 
lytically. In the star-connected network, we discover the 
phenomena of two-step phase transition. Through the 
simplification, we get the mechanism of this transition 
through the similar method as the case for a single star 
network, namely ensemble OA method. A linear rela¬ 
tion between initial order parameter and critical coupling 
strength between stars are get and checked. 


V. CONCLUSION 

In this paper, the star network and star connected net¬ 
work are used to study the main character of the het- 
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